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TECHNICAL MEMORANDUM 


DESIGN OF A MULTIPLE-PLY LAMINATED COMPOSITE TAPERED BEAMS 

INTRODUCTION 

In a recent paper, 1 the author developed a method based on classical lamination theory (CLT) 
and a linear elastic solution for the optimum design of a special class of laminated cantilever beams. 
The method, developed for beams which are tapered in width and thickness, expresses the exten- 
sional and bending stiffness matrices as functions of the linearly varying thickness along the length 
of the beam. The method also expresses the bending moment per unit width as a function of the 
linearly varying width along the length. The special case analyzed was the case of a four-ply 
laminated beam where the two outermost plies were tapered. The present work extends the original 
research to laminated beams of multiple plies. 

Findings of the present research include very high dependence of the interlaminar shearing 
stresses on the difference in layup angle between a tapering ply and the one adjacent to it. Other 
findings are that due to the strong dependence of the bending stresses on the stiffness parameter 
Du, the location of the highest bending, and consequently transverse (interlaminar) shear stresses, 
can be determined by knowing the shape of the Du curve as a function of beam length. It is also 
found that in some instances, stresses can be reduced by over 50 percent by simply interchanging 
plies within the design laminate with very small changes in the corresponding beam stiffness (less 
than 8 percent). The work presented here serves as a precursor to future research where the effects 
of the time-dependent behavior of the viscoelastic matrix of the composite material is included. As 
the use of composite materials becomes more widespread for space applications, such as for the 
Advanced X-Ray Astrophysics Facility-S (AXAF-S) project, design engineers will need more 
knowledge of long-term structural behavior of composite systems. 

LAMINATED BEAM EQUATIONS 

Following the CLT, the assembled extensional and bending stiffness matrices for an 
orthotropic material can be expressed as, 

-All A 12 Ai6 

[Aij] = A 2 1 A 22 A 2 6 (1) 

_ A^i Aft 2 A(56 _ 

Dll Di2 Di6 

[Dy] = D 2\ &22 D 2 6 (2) 

_ D61 D62 D66 _ 

The individual terms of equations (1) and (2) can be obtained following the procedure pre- 
scribed by the CLT with the exception that only the innermost plies are not tapered and their thick- 
ness does not depend on x. The outer plies are tapered and must be expressed as functions of x. In 
practice, plies cannot be tapered. In this report, the length of one tapered ply is equivalent to the 
length of several stepped plies which simulate the slope of the desired taper (figs. 1 and la). With 
this in mind, one can write. 
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Figure 1. Laminated tapered beam configuration. 





where Qf are the transformed reduced stiffnesses for the kth ply. Plies k m through m are the tapered 
plies at die bottom of the laminate. Plies k n through n are the tapered plies at the top of the laminate. 
L represents the length of the ply, and t k is the thickness of the ply. The validity of equations (3) and 
(4) is demonstrated in reference 1. One should notice that in equations (3) and (4), the terms 

[zrt(‘-r)] * nd hW'-r)] 

account for the dependence of the stiffness parameters Ay and Dy on the linear variation of thickness 
along the length of the beam. 


3 



For the cases studied, only pure bending of tapered beams is considered. With this assump- 
tion, the extensional stiffness matrix will not play a role in the solution of the beam problem. This 
would not be the case if the solution were to include the effects of transverse shear deformations. 2 
The constitutive relationships for the beam subjected to pure bending can be expressed as. 



^11 D i2 ^16 


D 


21 


Az2 D 26 


lA,, d, 


'62 ^66 


-1 


M x 

M y 

\M, 


(5) 


*y , 


For the beams studied, only bending along the longitudinal axis of the beam will be con- 
sidered. This implies that 


My = M xy = 0 


(6) 


The strain-displacement relationships for the typical Euler-Bemoulli beam are expressed as 


y,=- 


8 2 w 

dx 2 


r v = - 


d 2 w .. „ d 2 w 


dy 2 


• 7xy = - 2 


dxdv 


(7) 


For beams with large length-to-width ratios, the deflection w can be treated as being a func- 
tion of x only. Equations (5), (6), and (7) can be combined to express the Euler-Bemoulli beam 
equation as, 

( 8 ) 


y _ (£ W _ jj* 
fx~ ,2 ~ u U m x > 

where D x x is the first term of the inverted bending stiffness matrix. Equation (2) can be inverted 

numerically and a value of Dj j can be obtained at each location of x desired along the length of the 
beam. Figure 1 shows the beam dimensions and geometry. 

Once the values for D* j are obtained, they can be easily curve fitted with a higher order 
polynomial to obtain an equation as a function of length. This equation is of the form 


= X r>'-> 


1=1 


(9) 


In the present formulation, a compliance function of the form of equation (9) is calculated for 
each beam segment with tapered plies. Figures 2 through 4 show the comparison of actual versus 
curve fitted values of D* x for three different beam configurations. The material and geometry 
parameters used are shown in table 1. 


Table 1. Material and geometry parameters. 


Geometry 

Material 

w e - 4.50 in 

E u - 21.0E6 lb/in 2 

w 0 - 6.00 in 

E 22 - 1.7E6 lb/in 2 

t e - 0.50 in 

C 12 - 0.65E6 lb/in 2 

t 0 - 1 .00 in 

v 12 -0.31 

L- 29.25 in 

V 2) -0.017 
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Figure 4. D' n versus length for a 12-ply tapered laminated beam [90/45/30/15/0/0]*. 

In equation (8), M x is the moment per unit width of the beam. Since the width of the beams 
under study vary linearly as a function of length, this term becomes, 

M - M. - **«•-*> nm 

* w x x+w 2 W x x+W 2 ’ ( 1U 1 

where 


P = applied load at free end = qw e 
q = uniformly distributed applied line load at free end 
w e = beam width at free end 
w 0 = beam width at fixed end 


w, = 




w, = w n 


substituting equations (9) and (10) into equation (8) leads to 


dx 2 


= & r^c '- 1 

«=i r 


p(L-xy 

H',JC+W 2 


(Ha) 

(Hb) 


( 12 ) 
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Performing the required multiplications, equation (12) can be expressed in polynomial form. 
The resulting equation is easily integrated to obtain slope and deflection values. 


where 


d_w_ 

dx 2 


= % A,*' -1 


i= 1 


r,P(L-x) 

A I “ w l x+w 2 ' 


(13) 

(13a) 


Integrating equation (13), one can obtain values for slope and deflection at any point along 
the beam. Since the origin of the beams under analysis is located at the base (fixed end) of the beam, 
the boundary conditions of zero slope and deflection at the fixed end allow the elimination of con- 
stants of integration. Thus the slope and deflection of the beam can be expressed as 


dw _ y A ,X* 
dx ~ k\ i ’ 



A ,x M 
0 + 1 )* * 


(14) 

(15) 


Equation (13) represents a smooth continuous function for isotropic materials and for lami- 
nated materials where all ply angles are equal. However, because the bending stiffness D n is a 
function of material properties as well as laminate layup and thickness, marked discontinuities can 
occur at locations where plies “taper off’ to zero. At these locations, the curvature can be discon- 
tinuous because it is calculated independently for each beam segment with tapered plies. The physi- 
cal problem, however, does not allow discontinuities in slope at a location where a ply thickness 
becomes zero. With this knowledge, the degree of the polynomial in equation (13) can be selected 
such that the curvature over the length of the entire beam is reasonably approximated. As an 
example, an eight-ply laminated tapered beam [90/45/-45/0] Jt in which only the interior 
[-45/0/0/-45] plies remain of constant thickness, produces the curvature plot shown in figure 5. 
Material properties and geometry are shown in table 1. The marked discontinuity occurs because the 
outermost plies taper down in thickness to zero at this location (x = 14.625 in); and because, again, 
for each segment, the curvature is calculated independently. Also, at this location the laminate at x 
greater than 14.625 in only has six plies with the outermost plies having a layup angle of 45°. The 
laminate at x less than 14.625 in has eight plies with the outermost plies having a layup angle of 90°. 
Therefore, the beam stiffness is discontinuous at this location (x = 14.625 in). By approximating the 
curvature with a polynomial, a continuous function which can be easily integrated will be available. 
Figure 5 also shows the approximation from equation (13). Equation (15) is plotted in figure 6 for the 
[90/457-45/0]^ beam. This figure shows the calculated deflected shape of the cantilever beam. 


STRESS CALCULATIONS 

The stress-strain relationships can be expressed in the classical sense as, 

<*,>«= • < 16 > 

where the curvatures yj are calculated as in equation (13). Since the problems investigated are pure 
bending, the strains e; become negligible, and the stresses become principally functions of the curva- 
tures. Figure 7 shows the normalized bending stresses for a laminated beam where all plies are 
isotropic. 
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Figure 7. Normalized bending stress versus beam length (multiple isotropic layers). 

This serves the purpose of illustrating the stress variation of each ply as a function of length. 
For this example, the outermost three plies are involved in the tapering of die beam. One should 
notice how the inner ply stress equals the outer ply stress at the point where the outer ply tapers to 
zero thickness. In other words, the stresses for plies 1 and 2 are equal at x = 9.75 in; and the 
stresses for plies 2 and 3 are equal at x = 19.5 in. This is to be expected since the bending stiffness 
is a continuous function due to the isotropic nature of the plies. One should notice that ply 1 tapers 
off at x = 9.75 and ply 2 tapers off at x = 19.5 in. The stress for ply 4 is shown to illustrate the varia- 
tion in stress for a nontapering isotropic ply. 

For laminated orthotropic beams in which the bending stiffness is discontinuous at locations 
where plies taper down to zero thickness, the stresses will also have a marked discontinuity. In 
some instances, the inner plies may be subjected to higher stresses. Note that for isotropic beams, 
one would normally expect the highest bending stresses to occur at the outermost fiber. Figure 8 
shows the bending stress distribution for a [90/60/45/30/15/0], beam. Note that the outermost ply 
tapers down to zero thickness at * = 9.75 in. At this point, the second ply becomes the outermost ply 
up to x = 19.50 in. The third ply then becomes the outermost ply for the remainder of the length of the 
beam. 


It is also important to note that the 60° ply begins tapering at x = 9.75 and ends with zero 
thickness at x = 19.50. Likewise, the 45° ply begins tapering at x = 19.50 and ends with zero thick- 
ness at x = 29.25. For this laminated beam the maximum bending stress occurs at x = 19.50 on the 
0° ply above the center line of laminate symmetry (not shown in the figure). 
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The remaining in-plane stresses due to bending cr y and are also obtained by substituting 
equations (5) and (4) into equation (16) and setting £j = 0. In this manner, one obtains 

X+flg^i+fl8»y . U7) 

• as) 

Figure 9 shows the variation of o y for the top plies of the symmetric laminate investigated in 
figure 8. For an isotropic beam, these stresses would be negligible; but for the orthotropic laminate 
studied, the maximum G y stress is 13,050 lb/in 2 . This stress occurs at x = 9.75 on the 30° ply above 
the center line of laminate symmetry. One can also notice that the stress at x = 19.50 on the 60° ply 
has a significantly high compressive o y stress. This shows that buckling or compressive failure is 
possible at locations where one would normally expect to see only tensile stresses (isotropic case). 
Figure 10 shows the variation of for the top plies of the beam. Again, for the isotropic beam, one 
would expect these stresses to be negligible; but for this case, the maximum stress is 16,850 
lb/in 2 . This occurs at x = 9.75 on the 30° ply. 



Figure 8. Normalized bending stress versus beam length [90/60/45/30/ 15/0]*. 
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Ply 1 (tapered from 0.0 to 9.75 In) 

Ply 2 (tapered from 9.75 to 19.50 in) 
Ply 3 (tapered form 19.50 to 29.25 in) 
Ply 4 (nontapered) 


x-coordinate (Beam Length - in) 

Figure 9. Normalized in-plane stress versus beam length [90/60/45/30/15/0] 


Ply 1 (tapered from 0.0 to 9.75 In) 

Ply 2 (tapered from 9.75 to 1 9.50 In) 
Ply 3 (tapered from 19.50 to 29.25 In) 
Ply A (nontapered) 





x-coordinate (Beam Length - in) 

Figure 10. Normalized in-plane shear stress versus beam length [90/60/45/30/15/0] 
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INTERLAMINAR SHEARING STRESSES 


Laminated composite materials have relatively weak interlaminar shear strength, the load 
being sustained primarily by the resin or matrix at the ply interfaces. The CLT assumes that points 
lying initially on a normal-to-the-middle plane of the beam remain normal-to-the-middle surface of 
the beam after bending. This is one of the same assumptions used in the development of small 
deflection theory of thin plates. 3 This assumption is equivalent to ignoring shearing strains in planes 
perpendicular to the middle surface of the laminate. This also means that the entire laminate is 
treated as a single continuous layer with linear through-the-thickness strain variation. This is 
known as the Kirchoff hypothesis. 4 For laminated beam analysis, this implies that the interlaminar 
shear stress o xz cannot be calculated from the laminate stress-strain relationships. Instead, we 

must rely on the equations of motion (equilibrium) of classical theory of elasticity for the determina- 
tion of these stresses. 


The first equation of motion can be expressed in terms of Cartesian stress components as 


da, da„ da. 


Bx By dz 


+ pb x = pf x , 


(19) 


where b x are external body forces and f x are acceleration forces. For the static problems analyzed, 
these forces will be neglected. The resulting equation of equilibrium becomes 


da r dct da 

— + — ^ + -^=0 . 
dx By Bz 


( 20 ) 


Since the beams under investigation are loaded with an evenly distributed load at the tip and 
since the length-to-width ratio of these beams is large, the variation of the in-plane shearing stress 
across the width dimension y is negligible. Equation (20) can now be expressed as 


da da 
— + — = 0 . 
Bx Bz 


( 21 ) 


Equation (21) can be solved for o xz in the following manner: 


f 

°xz = ~ —dz. ( 22 ) 

Jz dx 

One should notice in equation (22) that if a x is a linear function of x, the derivative BaJBx is 
constant. This implies that a xz would be a function of the thickness coordinate z only . This would be 
the case if the cross sections of the beams under investigation were constant. Since the investiga- 
tion deals with beams of linearly varying width and thickness, the nonlinear variation of a xz as a 
function of length must be considered. Following Timoshenko 5 and Gere and Timoshenko, 6 the 
appropriate expression for g xz can be derived. Assuming that the beam of rectangular cross section 
in figure 1 1 does not have high rate of change in cross sectional thickness and width, the magnitude 
of the shearing stresses can be calculated by applying the same method used for prismatic beams. 
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Figure 11. Element of area for transverse shear stress determination. 





Summation of forces in the x direction leads to the following equation of equilibrium of the 
shaded element, 

(<•*<» i 

<7 x iw + dx = j 2 ( <J x +da x )(w+dw)dz - I 2 (c T x w)dz . (23) 


In equations (23) through (36), w is a width variable and not a deflection. Expanding equation 
(23) and neglecting products of small differential quantities, one obtains 

r -2**l , . x r _<&H - (Q 

<j xz wdx = w J 2 ^(7 x + -^ dx^z + J 2 (a^wWz-j 2 (<Tj ( w)dz • (24) 

The expression for bending stress for the fcth ply is obtained in a similar manner as equations 
(17) and (18) 

<7® = Z (*>M4ex, + e!?D,- 2+ fi<x) • (25> 

The quantity in the parenthesis is expressed in condensed notation as (QD x )^ k \ and equation (25) is 
rewritten as 


<rf _ z ^M x (QD x ) ik) . 

The quantity in the parenthesis of the first integral in equation (24) can be expressed as 




The integral can now be expressed as 


(26) 


(27) 


w 





d[M x (QD x ) (k) 

dx 



(28) 


Performing the necessary integration and neglecting products of small differential quantities, 
one obtains the following expression: 


TT M x (QD x ) (k) + [M x (QD x ) (ky ] - w[ ^ >]2 [M x (QD x ) (ky ] + d W x (QD x ) (k} ] ^ 


d[M x (QD x ) <ky ] 
2 dx 


dx 


(29) 


dx 


Substituting equation (26) into the second integral on the right-hand side of equation (24) 
and performing the necessary integration yields the following expression after neglecting products of 
small differential quantities: 


-jp [M x (QD x ) iky ] - l z - (k) ^ dw [M x (QD x ) (ky ] . 


(30) 
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Following this same procedure, the third integral on the right-hand side of equation (24) 
becomes 

[M X (QD X ) (k) ] - [z( ^ 2w - [M x (QD x ) (ky ] . (31) 


ing: 


to 


Adding expressions (29), (30), and (31) according to equation (24), one obtains the follow- 


<? xz wdx = -j- [M x (QD x y >]dt + -g ^ dx 5 dx 

+ [M x (QD/ ky ]dw - [M x {QD x ) (ky ]dw . 

Solving for cr^, the following is obtained: 


(32) 




r r M(nr >^ * Jt 2 [z<*>3 2 V ^(gi>,) (t> l , MJQPJ™ dw] 
j[M x (QD x ) % + ^-g 5 r\ dk + w ^ 


(33) 


For the beams under analysis, both the width and thickness are linear functions of x according 


t(x) = t i x+t 2 , 

w(x) = w x x+w 2 , 


(34a) 

(34b) 


where 


‘1 “ / 


W, = 


L 

W e -w b 


h = h . 

w 2 = w b . 


T alcin g the derivatives of equations (34) with respect to x leads to 

dt _ , _ te-tb 

dk- 1 '- L ’ 


dw _ 


(35a) 

(35b) 


Substituting the results of equations (35) into equation (33) yields the final expression for 
interlaminar shear stress for the kth ply. 


^ = h^lA [Mxi QD/ k >] + 


0 1 *** 2 ) 


8 


--*?] { 


d[M x {QD x ) (ky ] , w, 


<£t 


WjX+Wj 


• (36) 


One should notice that equation (36) gives the value of o xz at any location along the length x 
and through the thickness z of the laminated beam. Due to limitations in the validity of equation (36) 
for large values of width-to-thickness ratio 7 8 9 and because for small values of this ratio, stress 
singularities at the free edges can cause delaminations; 10 11 our studies are limited to 
6 < w/t < 10. 
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As observed by Whitney, 11 the interlaminar shear stress variation through the thickness of 
the laminate is highly dependent on the lamina stacking sequence. Figure 12 shows a comparison 
between the transverse shear stress for two orthotropic beams; [90/45/30/ 15/0/0], and [0 6 ],. Both 
beams have the same geometry and material; however, only the [0 6 ], laminate produces the stress 
distribution which is typical for the fixed end of a doubly tapered cantilever beam. If the material were 
isotropic, this stress distribution would be similar to the [0 6 ], laminate. Figures 13 through 15 show 
these stresses as calculated at various locations along the length of the beam. One should notice 
how, in all cases, the [90/45/30/15/0/0], develops considerably higher stresses. It is also of interest 
to note that at x - 29.25, where the tapered beam consists of only [0 2 ], plies, the shape of the 
distribution follows the classical sense for both laminates. This is to be expected since there is no 
stiffness mismatch between plies at this location. 



Figure 12. Transverse shear stress comparison for tapered beams 
(orthotropic versus isotropic). 
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z-coordinate (Beam Thickness) at x s 9.75 in 

Figure 13. Transverse shear stress comparison for tapered beams 
(orthotropic versus isotropic). 



Figure 14. Transverse shear stress comparison for tapered beams 
(orthotropic versus isotropic). 
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z-coordinate (Beam Thickness) at x = 29.25 in 

Figure 15. Transverse shear stress comparison for tapered beams 
(orthotropic versus isotropic). 


As a direct result of the variation in stiffness along the length of the beam and because 
certain plies taper off to zero thickness at specific locations, there can be significant differences in 
stress on either side of the point where a ply “tapers Off.” It has been seen that the greater the 
difference between the ply layup angle of two adjacent plies, the greater this difference in stress. 
Figure 16 shows a comparison of the outermost ply bending stresses for three laminates. These are, 

a. [45 5 ] s 

b. [45/90/80/60/45] s 

c. [ 80/90/45/60/45] s . 

The geometry and material properties are the same as the beams previously investigated. 

One should notice that the only difference in b and c is the interchange of the 80° and 45° plies. 

These plies have a difference in ply layup angle with the 90° ply of 10° and 45° respectively. It is evi- 
dent that for both cases b and c, the greatest difference in stresses at the point where a ply tapers off 
occurs at the 45/90 or 90/45 interface. Also, the magnitude of the largest value of stress is consider- 
ably smaller for the c laminate. This is an important observation since the overall spring rate for both 
the b and c laminates changes very little. A summary is presented in table 2. 
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Figure 16. Outermost ply bending stress comparison for three orthotropic laminates. 

Table 2. Summary of stress and stiffness for tapered beam 
(top ply only) (200 lb applied load). 


Laminate 

Max. Bending Stress 
(lb/in 2 ) 

Spring Rate 
(lb/in) 

[ 45 s ], 

6,821 

53.62 

[45/90/80/60/45], 

19,711 

56.98 

[80/90/45/60/45], 

8,010 

61.35 


It is seen that a reduction in maximum bending stress of 59 percent occurs by simply changing 
the location of the 45° and 80° plies. The spring rate, however, increases by only 7.6 percent. Figure 
17 shows a comparison plot of the bending stiffness component D n for the b and c laminates. It is 
relatively straight forward to see that the effect of the tapering plies on the b laminate is consider- 
ably greater. This is evident by the higher rate of change of the stiffness over the first 10 inches of 
the beam length. Figure 17a shows the bending stiffness components D n , D n „ and D l6 for laminate 
b. The same discontinuity at x = 9.75 in is evident for all three curves. 
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Figure 17a. Comparison of bending stiffness components D n , D n , and D 16 [45/90/80/60/45],. 
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The D,y stiffness plots can be used as an indicator for determining at which locations stresses 
will be highest. The following rules are proposed: 

Stiffness/Stress Rules 


• £>i i curves with highest rate of change indicate 

highest stresses 

• Sharp discontinuities, “knees,” in the D n curves 

indicate locations of highest stresses. 


As an illustration, figure 18 shows the bending stress distribution along the length of the 
beam. One can easily notice that at the point where the top ply (45° ply) tapers off, a sharp discon- 
tinuity in stress occurs. A corresponding stress change occurs at the 90°, 80°, and 60° plies at this 
location. 

The change in stiffness at the point where the second ply tapers off is almost negligible; 
likewise, the change in stress at this location ( x = 19.50 in) is small (smooth transition). 

From equation (22), it is evident that the transverse shear stress a xz is a direct function of 
the bending stress a x . With this knowledge, using the proposed stiffness/stress rules, one can 

analogously predict where the maximum interlaminar shear stresses will occur. Figure 19 shows 
these stresses for the three outermost (tapered) plies of the beam. One should notice how the 
maximum stress occurs, as expected, at the point where the first ply tapers down to zero thickness. 
This is also the location where the greatest stress discontinuities occur. These discontinuities, 
again, are due to two reasons: 

• The outermost plies are tapering which contribute to a rapid increase in the compliance 
matrix (decrease in stiffness matrix). At the point where the ply tapers, the stiffness 
change is analogous to that of an isotropic beam with a sudden change in cross-sectional 
area. 

• The difference in the layup angle between the tapering ply and the adjacent ply is 
greatest. 

Another parameter which provides meaningful information about the relative magnitude and 
location of highest stresses is the nondimensional stiffness parameter ( QD X ) [see equations (25) 
and (26)]. Figures 20 through 23 show how this parameter varies as a function of length of the beam. 
Figure 20 clearly shows a marked difference between the first (45°) and second (90°) plies while 
there is a small change between the second (90°) and third (80°) plies. In this figure, all plies are 
tapered which means they are the outermost plies. In figures 21 and 22, one can notice the large dif- 
ferences in ( QD X ) for the 80-to-60 and 60-to-45 ply interfaces respectively. Large stress discon- 
tinuities would be expected at these plies based on this notion. However, one should recall that in 
figure 17, the bending stiffness parameter D n had the highest rate of change up to the x = 9.75 in 
location. Between 9.75 and 29.25 in, the parameter had a very smooth curve. This implies that in 
order to predict the relative location and plies that have the largest stresses, it is important to know 
the difference in the ( QD X ) parameter between adjacent plies and the rate of change of the D n bend- 
ing stiffness parameter as a function of beam length. With this in mind, the magnitudes and discon- 
tinuities of the interlaminar stresses in figure 19 can be explained. 
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x-coordinate (Beam Length - in) 

Figure 18. Bending stress distribution along the length of a tapered orthotropic laminated beam 

[45/90/80/60/45],. 



Figure 19. Interlaminar shear stress variation for tapering plies [45/90/80/60/45],. 
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Figure 20. [QD x ] k parameter for laminated beam with tapered orthotropic plies 

[45/90/80/60/45],. 



Figure 21. [ QD x ] k parameter for laminated beam with tapered orthotropic plies 

[45/90/80/60/45],. 
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x-coordinate (Beam Length - in) 

Figure 22. [ QD x ] k parameter for laminated beam with tapered orthotropic plies 

[45/90/80/60/45],. 
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Figure 23. [QD x ] k parameter for laminated beam with tapered orthotropic plies 

[45/90/80/60/45],. 






Since the bending stiffness parameter D n has the highest discontinuity at x = 9.75 in, this 
should be the location of the highest overall stresses. This can clearly be seen in figure 19 where the 
differences in <7 XZ at the x = 19.50 in location are considerably lower. 

Figures 24 and 27 are plots of the interlaminar shear stresses between adjacent plies at the 
x = 9.75 in and x = 19.50 in location. These plots also show the ply angle corresponding to each 
lamina. Figure 24 shows a large stress difference between the 45° ply and the 90° ply at the 
x = 9.75 in location while the stress difference between the 90° and 80° plies is almost negligible. 
This corresponds directly to the difference in the ( QD X ) parameter in figure 20 at the same locations. 
Figure 25 indicates a larger difference between stresses at the 80° to 60° ply interface at 
x = 19.50 in. Again this compares directly with the ( QD X ) parameter difference in figure 21. The 
magnitude of these stresses, however, is lower because the £>j*j curve has approximately the same 
rate of change on either side of the x = 19.50 in location. The difference in stresses in figure 26 again 
corresponds to the difference in ( QD X ) in figure 22. One should notice, however, that the stress 
between the 60° and 45° ply does not go to zero at the free end (x = 29.25 in) of the beam. This is 
because at the free end, where the load is applied, the transverse shear stress distribution allows 
stresses to be zero only at the outermost plies (80° to 60° interface) and not necessarily at any of 
the inner plies. Again figures 27 and 23 show the same relationship between a xz and ( QD X ) pre- 
sented in the previous plots. 


SUMMARY 


The design engineer must have knowledge of the effects of material and geometric properties 
on the magnitude and location of stresses and deformations on designed hardware during the design 
process. This is particularly true when designing complex components fabricated from composite 
materials. For the case of laminated tapered beams, the following recommendations are offered: 

• The difference in ply angle between a tapering ply and the adjacent one should be as small 
as possible to reduce the magnitude of the bending and interlaminar shearing stresses and 
their discontinuities at the point where the ply “tapers off.” 

• A laminate should be designed so that it will produce the smallest rate of change in the 
bending stiffness components D The laminate design should produce the smoothest, 
although not necessarily, continuous curve possible. 

• The designed laminate should exhibit the smallest differences in ( QD X ) between tapering 
and nontapering plies. 

In reading these recommendations, it is obvious that a completely isotropic beam will meet all 
of them. It is also noted that laminates in which all ply angles are the same also meet these recom- 
mendations. Figure 28 shows a plot of the bending components D n , D 12 , and D l6 , for a [30 5 ] s 
laminate. There are no discontinuities in any of the curves indicating a smooth transition of stresses 
between tapering and adjacent plies. One should notice however that the rate of change of these 
curves decreases towards the tip of the beam. This along with the fact that the bending moment per 
unit width (fig. 29) increases toward the fixed end of the beam indicates that the highest stresses 
are to be expected toward this end. Another interesting plot is figure 30 where ( QD X ) plots for the 
[30 5 ] s laminate are all identical. This is to be expected since all plies have the same transformed 
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reduced stiffnesses and the stiffness terms do not produce discontinuous curves. Figure 31 
shows the bending stress distribution for this laminate. It is important to note that it has the same 
kind of stress distribution exhibited by an isotropic multiple ply laminate (fig. 7). 

Efficient and effective structural components can be manufactured with the knowledge of 
material and physical characteristics of the design. This study offers these design considerations for 
a particular class of laminated beams. Future work could possibly involve beams of more complex 
cross sections such as those typically used in structural design applications. Cross sections such as 
I, C, or H shapes can twist when bending loads are applied. Tapering such beams may result in 
increased warping resistance and decreased interlaminar shearing stresses. Structural optimization 
of the laminate can also minimize the weight while making the member structurally sound. 
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